Gap size and capture zone distributions in one-dimensional point island nucleation and growth simulations 



o 

(N 
o 

Q 



O 

< 



> 

cn 
m 



X 



Gap size and capture zone distributions in one-dimensional point island 
nucleation and growth simulations: asymptotics and models 

K. P. 0'Neill/'E3 M. Grinfeld,i'0 W. Lamb/'El and P. A. Mulheran^H 

Department of Mathematics and Statistics, University of Strathclyde, Glasgow 
^''Department of Chemical and Process Engineering, University of Strathclyde, 
Glasgow 

(Dated: 12 December 2011) 

The nucleation and growth of point islands during submonolayer deposition on a one-dimensional substrate is 
simulated for critical island size z = 0, 1, 2, 3. The small and large size asymptotics for the gap size and capture 
zone distributions (GSD and CZD) are studied. Comparisons to theoretical predictions from fragmentation 
equation analyses are made, along with those from the recently proposed Generalised Wigner Surmise (GWS). 
We find that the simulation data can be fully understood in the framework provided by the fragmentation 
equations, whilst highlighting the theoretical areas that require further development. The GWS works well 
for the small-size CZD behaviour, but completely fails to describe the large-size CZD asymptotics of the 
one-dimensional system. 
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I. INTRODUCTION 

The nucleation and growth of islands during submono- 
layer deposition is of considerable theoretical interest as 
a fundamental problem in the statistical mechanics of 
growth processesi"— . The sizes and spatial organisation 
of the nucleated islands ultimately determine the higher- 
level structures such as film and nanostructure array 
morphologies^. A long-established strategy in the analy- 
sis of the statistical properties is to study capture zones 
of islands, since these not only reflect spatial organisa- 
tion but also determine growth rates of the islands^r— . 
Therefore, the evolution of capture zones during the de- 
position process has been a focus of many recent theo- 
retical works^"— . 

Recently, Pimpinelli and Einstein introduced a new 
theory for the capture zone distribution (CZD) employ- 
ing the Generalised Wigner Surmise (GWS) from random 
matrix theoryi^, causing some controversy. Oliveira and 
ReisJ-7 have presented simulation results for islands grown 
on a two-dimensional substrate with critical island size 
i — 1 and 2, providing some support for the proposed 
Gaussian tail of the CZDi^. However, Li et. al}'^ pre- 
sented an alternative theory which yields a modified form 
for the large-size CZD behaviour, supported by data for 
the simulated growth of compact islands with i = 1. This 
form seems to agree with that found by Oliveira and Reis, 
contradicting the GWS"'^^. In other work, Shi et. alJ^ 
studied i = 1 models in o? = 1,2,3,4 dimensions, find- 
ing that the CZD is more sharply peaked and narrower 
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than the GWS suggests. Therefore it is by no means es- 
tablished whether the GWS provides a good theoretical 
basis for understanding the distribution of capture zones 
found in island nucleation and growth simulations. 

The simplified case of point island nucleation and 
growth in one dimension has proven to be a good test 
case for theories. For example, Tokar and DreysseiS have 
recently used this model to illustrate their accelerated ki- 
netic Monte Carlo algorithm for diffusion limited kinet- 
ics, finding excellent scale invariance in the island size 
distribution. Blackman and Mulheran^^ studied the sys- 
tem with critical island size i — 1, using a fragmentation 
equation approach. In this system, we can view the sub- 
strate as a string of inter-island gaps, and new island 
nucleation caused by the deposited monomers as a frag- 
mentation of these gaps. Thus in order to understand 
the CZD, it is important first to be able to describe the 
gap size distribution (GSD). 

In recent work^S we have extended the analysis of the 
fragmentation equations of [l3 to the case of general i = 
0, 1, . . .. We have been able to derive the small and large 
size asymptotics of the GSD, and by assuming random 
mixing of the gaps caused by the nucleation process, we 
have also derived the small size asymptotics for the CZD 
for general i and the large size behaviour for i = 0. 

One key feature to emerge from this work is that the 
asymptotic behaviour of the CZD is again different to 
that of the GWS^'^. It therefore is appropriate to ask 
what support, further to that in reference [l^, for the 
fragmentation equation approach is offered by Monte 
Carlo simulations of the system. Recent work by Gonza- 
lez et. al^ has revisited the case of i = 1, developing the 
original fragmentation equatioiii^ and GWS arguments 
in response to deviations between prediction and simula- 
tion. In this work we will explore simulation results for 
the one-dimensional (1-D) model with i — 0,1,2,3, and 
consider the relative merits of the fragmentation theory^S 
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and GWS-"^ approaches. 

The paper is organised as follows. In Section II we 
summarise the relevant theoretical resultsi^i^S. In Sec- 
tion III we describe the Monte Carlo methods used in 
our work, both for the full simulation of the island nu- 
cleation and growth processes as well as for nucleation 
within single gaps. Simulation results are presented in 
Section IV and compared to theoretical predictions, and 
we finish with a summary and our conclusions in Section 
V. 



II. THEORY AND PREDICTIONS 
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FIG. 1. Summary of the features of the one- 
dimensional point island model with i — 1. 
Solid circles represent an island; open circles 
are monomers. A capture zone is the separa- 
tion of the bisectors of neighbouring gaps.— 



The data from MC simulations can be used as a bench- 
mark against which to test predictions of theories for the 
GSD and CZD. In the next two subsections, we will dis- 
cuss the predictions of two competing theories, namely 
the fragmentation equation approach and the GWS. 

A. Fragmentation Equations 

We follow the Blackman and Mulheran approach for 
the 1-D point-island model with i — I-IS. Island nu- 
cleation events are viewed as the fragmentation of gaps 
between stable islands; see Figure [T] A nucleation that 
occurs in a parent gap of width y will result in the cre- 
ation of two daughter gaps of widths x and y — x. The 
probability that the nucleation occurs at position a; < y is 
taken from the long-time (steady state) monomer density 
profile in the gap, 

ni{x) = ^x{y - x), (1) 

where R = D/F is the ratio between the monomer dif- 
fusion constant D and the monomer deposition rate F . 
In particular, we assume that the nucleation probabil- 
ity is obtained from this monomer density ni{x)°' , with 
the value of a reflecting the nucleation process. We then 
obtair*2£ 

d 

^u{x, t) = -B(a + 1, a + l)x^u(x, t) 
at 

+ 2 [x{y-x)ru{y,t)dy, (2) 

J X 

where B{-,-) is the Beta function and A = 2a -I- 1. Here, 
u{x,t) is the number of gaps of size x at time t. The 
first term on the right hand side of ([2]) is the rate at 
which gaps of size x are removed from the population 
by a nucleation event. The second term describes the 
creation of gaps of size x from the fragmentation of larger 
gaps, with the factor 2 reflecting the symmetry of the 
fragmentation kernel. 

In [23] we set a = i + 1 under the assumption that 
nucleation is a rare event solely driven by the diffusion 



of the monomers. In doing this we implicitly assume 
that the i + 1 monomers necessary to create the nucleus 
are all in some sense mature, each separately obeying the 
long-time steady-state density profile ni{x). However, we 
shall also have need to consider the case when nucleation 
is triggered by a deposition event. Here a newly deposited 
monomer either lands close to (or even directly onto) a 
pre-existing cluster of i > mature monomers. In this 
case, we set a = i. 

Eqn. ([2|) admits similarity solutions of the ioimM^^— 

u{x,t) = x{ty^(t>{x/x{t)), (3) 

where x{t) is the average gap size. The following asymp- 
totics are then founci^^: 

(/)(z) ~ fcz" as z ^ 0; (4) 

(/)(z) ^ exp(— cz'*') as z ^ oo, (5) 

for constants c > and k. Here z — x / x{t) is the scaled 
gap size. 

We may use this information to understand the scal- 
ing asymptotics of the CZD P(s) where s is the scaled 
capture zone size. On a 1-D substrate, a point island's 
capture zone is made up from half of the gap to its left 
combined with half of its gap to the right (see Fig. [TJ . 
If there is no correlation between the sizes of two neigh- 
bouring gaps, we can writei^ 



P{s) = 2 / 0(z)0(2s - z)dz. (6) 

^0 

The factor 2 is included to preserve the normalisation for 
P(s). The small-size asymptotics of P{s) is thenM 

P{s) - fcs^^+i as s ^ 0, (7) 

for some constant k. The large size scaling of P{s) can 
be computed only for the special case a = l,i = 0. It 
has been shown that-", for some constant fc. 
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P(s) - fcs-^/2g-2sVA.^ as s oo, (8) 

where /i is a positive constant. 

We note here that the large-size asymptotics of the 
GSD and the CZD are thus the same for spontaneous 
nucleation. Given the form of Eqn. © for P{s), we con- 
jecture that the correspondence between the GSD and 
CZD large-size asymptotics will hold for other values of 
a — 2,3, 4... although it has not been proved. 

B. Generalised Wigner Surmise 

Recently, Pimpinelli and Einsteini^ conjectured that 
the CZD is well described by the Generalised Wigner 
Surmise (GWS) formula, which depends only on one pa- 
rameter P that reflects the critical island size i and the 
dimensionality d of the substrate: 

Ppis)^aps^expi-bps''), (9) 
where /? is given by 

Here and 6^ are normalisation constants so that 

/ Pp{s)ds = / sP0{s)ds = 1. 
^0 "'0 

The remarkable feature of this conjecture is its univer- 
sal nature; unlike the fragmentation equation approach 
described above, which is specific to the 1-D substrate, 
the GWS is claimed to hold for all dimensions. Pimpinelli 
and Einstein^ demonstrated good agreement with simu- 
lation results taken from the litcratur o^°i^^ , but only with 
j = 1 in d = 1. In their most recent work^i, this group 
analyse the i=l, d=l model in more detail and modify 
equations ([9]) and ([T0| in response to their findings. Here 
we note that the asymptotics of the fragmentation equa- 
tion approach above and the GWS do not agree-" for all 
i, whether we adopt a = iora = i + l. This in part 
motivates the present Monte Carlo study. 

III. MONTE CARLO SIMULATIONS 

We perform Monte Carlo simulations of point islands 
on a 1-D substrate. We adopt the same methodology 
as in previous work for i = 1—, but here we now also 
simulate a range of values of the critical island size 
i = 0,1,2,3. In the first subsection we describe the 
full simulation for island nucleation and growth, and in 
the second we describe a variant for obtaining nucleation 
rates within a single gap. 



A. FULL SIMULATION 

In the full simulation, monomers are randomly de- 
posited at rate F monolayers per unit time onto an ini- 
tially empty one-dimensional array of sites representing 
the substrate. Deposited monomers diffuse at the rate D 
by performing random hops between nearest neighbour 
lattice sites. We use periodic boundary conditions. For 
i > 0, if the monomer number at any site exceeds the crit- 
ical island size, a new island is nucleated. In the case of 
spontaneous nucleation {i = 0), monomers have a small 
probability p„ of nucleating a new island each time they 
hop. Once nucleated, an island increases in size by ab- 
sorbing any monomer which hops onto it from a nearest 
neighbour site. In the work discussed here, the islands 
only ever occupy one lattice site whatever their size in 
absorbed monomers. These processes are illustrated in 
Fig. [II 

As the deposition rate F increases, the average time a 
monomer diffuses before meeting another monomer de- 
creases. Due to the competition between diffusion and 
deposition, the statistical properties depend on the ratio 
R = D/F. 

The nominal substrate coverage, 6 = Ft, is a useful 
measure of the extent of the deposition process. Note 
that because we simulate point islands, this coverage can 
be greater than 100% even whilst most of the substrate 
remains free for monomer diffusion. For a fixed value 
of 9, the average distance between islands increases if R 
is increased. Similarly, for fixed R, as coverage increases 
the island density also increases. We are interested in the 
scaling properties of the aggregation regime^, where the 
island density exceeds the monomer density. The value 
of 9 for which this regime starts depends on i and R, and 
we check that the values of 9 employed are sufficiently 
high to ensure that we are in the aggregation regime. 

Our simulations were performed on lattices with 10^ 
sites, with R = 8 x 10^ up to coverage 9 = 100%, aver- 
aging results over 100 runs. For ? = we set the sponta- 
neous nucleation probability to p„ = 10~^. With these 
parameters, we find island densities of about 0.5%, 1.5%, 
0.5% and 0.25% for i = 0, 1, 2, 3 respectively at 61 = 100%. 
We note that this is a long way short of the limit referred 
to by Ratsch et. al?^ where scaling breaks down as the 
lattice becomes saturated with islands. We also have no 
finite size effects with this size of lattice, and do not need 
to implement accelerated algorithms'^. 

B. SINGLE-GAP NUCLEATION RATE SIMULATION 

In the single-gap simulation, we simulate island nucle- 
ation events in gaps ranging from size g — to g — 500, 
which proves to be adequate to illustrate the nucleation 
mechanisms at play. In this variant, monomers can dif- 
fuse as usual on a lattice of length g, but are removed 
from the simulation if they try to hop beyond the ends 
of the lattice. 
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We set the nominal monolayer deposition rate F to 
unity, so that a monomer deposition increments the sim- 
ulated time by 1/5 (recall that it is the ratio R — D/F 
which is important, rather than the absolute value of 
either F or D). Upon each step of the algorithm, we 
either deposit a new monomer at a randomly chosen 
site in the gap, or diffuse an existing monomer accord- 
ing to the relative rates of these processes. Explicitly, 
a monomer is deposited into the gap with probability 
Fxg/{Fxg + Dxn) = l/(l-|-i?x n/g), where n is 
the number of monomers currently in the gap. If no de- 
position occurs, a randomly chosen monomer hops to a 
nearest neighbour site. If i -I- 1 monomers (for i — 1,2,3) 
coincide at a site to form a stable nucleus, the simula- 
tion ends and the time to the nucleation event recorded. 
Repeat runs always start with an empty lattice, and are 
used to obtain reliable statistics on the nucleation times 
within each gap size. 

We use R = 10^ for i = 1, 2, and R = 10^ for i = 3 (due 
to simulation time constraints). Finally, we also moni- 
tor the average monomer density profile across the gaps, 
along with the number of hops each monomer makes 
in the single-gap simulations. The latter will indicate 
whether or not island nucleation is influenced by depo- 
sition events. If nucleation is caused solely by the dif- 
fusional fluctuations of monomers, then the stable nu- 
clei should only include monomers that have taken many 
hops. If however nucleation closely follows a deposition 
event, then the nuclei will contain monomers that have 
only made few hops since their deposition. 



IV. RESULTS 



large gaps can also be seen in our full simulation results 
(data not shown). This will have direct consequences for 
how the nucleation rate varies with gap size for larger 
gaps, as we now show. 
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A. SINGLE GAP NUCLEATION RATE 

In Fig. [5] we show the results for the average monomer 
density profile within gaps of size 5 = 100 and g = 300 
for i = 1. For the smaller gap size, we see that the profile 
agrees well with the assumption made in the fragmenta- 
tion equation approach^^, coinciding with the long-time 
steady-state solution of the diffusion equation with ran- 
dom deposition (Eqn. This is typical for the lower 
end of the range of gap sizes that occur in the full simu- 
lation at higher coverage, for all the values of i that we 
have studied. 

However, for the larger gap size g — 300 shown in 
Fig. [21 we see that the monomer density profile falls a 
long way below the long-time prediction. This behaviour 
is typical for all values of i at the upper end of gap sizes 
found in our full simulations. The reason for the shortfall 
is the higher nucleation rate in the larger gaps; the aver- 
age monomer density profile does not have sufficient time 
to reach its saturated level in Eqn. ([1]) before a nucleation 
event occurs. As stated, the range of gap sizes g used in 
the single-gap simulation is determined by the range typ- 
ically seen in our full simulations. Therefore, this failure 
to reach the saturated monomer density profile with the 
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FIG. 2. Monomer density profile in a single gap of size g — 100 
(top) and g — 300 (bottom) for i = 1 

In Fig. |3] we show the average time for a nucleation 
event to occur {tnuc) for all single gaps in the case of 
i — 1,2 and 3 (note that the data for i = 2 and i = 3 have 
been shifted horizontally to avoid overlapping curves). 
We note that the data obeys the power-law form pre- 
dicted by the fragmentation equation approach for small 
gap sizes g, but as expected deviates strongly for larger 
gaps. In fact, the average time to nucleation becomes 
much higher than predicted by the use of the saturated 
monomer density profile, since the actual profile for the 
larger gaps is lower, therefore presenting slower than ex- 
pected nucleation rates (but still fast compared to the 
time it takes for the monomer density to grow from zero 
to its saturation level). 

The straight line fits in Fig. [3] are for the small gap 
size data only {g € [50,150]). We use these to estimate 
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how the nucleation rate varies with gap size g through 
^/{tnuc) Of 9'^, with the values of the power 7 reported 
in Table m We have used bootstrap methods with 1000 
samples of size as big as 80% of the original to find an 
approximate 95% confidence interval in Table HI 

The fragmentation equation approach (Section IIA 
above) suggests that this power should be 2i + 1 or 
2i + 3, depending on whether island nucleation is driven 
by monomer deposition or solely by monomer diffusion. 
The results in Table |T] suggest that the simulation re- 
flects both these mechanisms, with the small gap size 
nucleation rate exponent lying between these two possi- 
bilities. 
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monomer density, helping to explain the intermediate val- 
ues for the nucleation rate exponents in Table HI 
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FIG. 4. Histogram of the number of hops taken by the 
youngest monomer in a nucleus for i — 1, for gap size g = 100 
(crosses) and g = 300 (diamonds). In the main figure the 
number of monomers is truncated at 100. The inset shows 
the same result at the lower number of hops without trunca- 
tion of the number of monomers in the histogram. 



FIG. 3. Average time for a nucleation event to occur at all 
gaps. 
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A" 


Simulation 


1 


3 5 


3.628 ± 0.035 


2 


5 7 


6.606 ±0.118 


3 


7 9 


8.283 ± 0.266 



A = 2j -I- 1 
A = 2i -f 3 

TABLE I. Small gap nucleation rate exponents from the single 
gap simulations. 



In Fig. 0] we present histograms for the number of 
hops taken by the youngest monomer in a nucleus for 
the g = 100 and g = 300 i = 1 simulations. The his- 
togram has a long tail, showing that in many cases all the 
monomers in the nucleus are indeed mature in the sense 
that they have diffused many times since their deposition. 
However, there is also a sharp increase in likelihood of a 
monomer only taking very few diffusive steps before be- 
ing caught up in a nucleation event. In other words, there 
are a significant number of nucleation events driven by 
fluctuations due to deposition. This supports the con- 
clusion that nucleation in these simulations is driven by 
a combination of deposition and diffusion fluctuations in 



B. FULL SIMULATION BEHAVIOUR 

Having established the nucleation behaviour in single 
gaps, we can now look at the results observed in our 
full Monte Carlo simulations. The fragmentation equa- 
tion approach again provides concrete predictions for the 
small and large size behaviours for the GSD and CZD. 
We will also be able to compare the CZD properties with 
the GWS, and establish which of the two theories pro- 
vides the better framework to understand the behaviour 
observed. 



1. SMALL SIZE SCALING OF THE GSD AND CZD 

In Figures [5] and [6l we report the small size behaviour 
of the GSD {(f){z)) and CZD (P(s)) in logarithmic scale 
a.t 9 = 20%. In order to fit the slopes in these plots, and 
obtain reliable error estimates, we adopt the following 
numerical technique. The size data are binned using reg- 
ularly spaced bins on the logarithmic abscissa, with bin 
widths 6™c where b and c are flxed constants and m > 0. 
By choosing a range of values for b — 1.1, 1.2, 1.3 and 
1.4, and c = 0.0125, 0.025 and 0.05, all of which pro- 
vide reasonable choices for binning the data, we obtain a 
number of straight-line fits. This allows us to calculate 
the average of these gradients and a 95% confldence in- 
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terval. The results of this fitting procedure are shown in 
Tables E] and imi 
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FIG. 5. Small-size GSD in logarithmic scale for i = 0, 1, 2 
and 3 at coverage 9 = 20%. The dashed line is the straight 
line fit to data. 
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2.719 ± 0.082 


3.271 ±0.056 



b a = i + l 
9 = 20% 

e = 100% 

TABLE II. Average gradient for the small size scaling of the 
GSD using different bin-widths at coverage 9 = 20% and 100% 



i 


2a ± r 2a ± 1^ 


GWS^ 


CZD^ CZD^ 





3 


2 


2.730 ± 0.030 2.751 ± 0.086 


1 


3 5 


4 


4.187 ±0.050 4.372 ±0.149 
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7.200 ±0.382 6.138 ±0.124 
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FIG. 6. Small-size CZD in logarithmic scale for i = 0, 1, 2 
and 3 at coverage 9 — 20%. The dashed line is the straight 
line fit to data. 

For the small-size asymptotic behaviour of the GSD 
and CZD we compare the data from MC simulations with 
the fragmentation equation approach predictions of Sec- 
tion IIA. For the GSD, the dominant term is as z ^> 
(see Eqn. (U)). Likewise, for the CZD the dominant term 
is s^"+^ (see Eqn. (O). For the latter, we also have the 
competing prediction of the GWS which is (Eqn. ([9])). 
The values from these theories are also displayed in Ta- 
bles HI] and HIl 

The results for the small size scaling exponent of the 



TABLE III. Average gradient for the small size scaling of 
the CZD using different bin-widths at coverage 9 — 20% and 
100% 



GSD in Table HI] show that the fragmentation equation 
approach provides a reasonably sound framework for un- 
derstanding the island nucleation and growth process. 
For i = 1, 2, 3 we see that the exponent at 9 = 100% lies 
between the two possible values a = i and a — i + 1 
suggested by the theory. This is as expected following 
the single-gap nucleation results presented above, which 
show that both the deposition- and diffusion-driven nu- 
cleation mechanisms are at play in the simulations. We 
note that the 9 = 20% results for i = 3 lie below 
a — i — 3, but we believe that this is due to the fact 
that the simulation has only just entered the aggregation 
regime in this case. We also see that for i = 0, the ex- 
ponent is close to the a ~ i + I — 1 prediction (a = is 
not a viable possibility), being closer at 6' = 100%. 

The trends shown in the small size scaling exponent 
of the CZD in Table Hill are rather similar. We see the 
i = data are close to the A = 2i ± 3 = 3 prediction of 
the fragmentation equation approach, being somewhat 
larger than the ;3 = 2(i ± 1) = 2 predicted by the GWS. 
For i = 1,2 the data are bracketed by the two alterna- 
tives suggested by the fragmentation theory, as indeed 
is the GWS exponent which appears to present a rea- 
sonable compromise given the two alternative nucleation 
mechanisms. The case of i = 3 provides an exception, 
which hints at the breakdown of the relation in Eqn. (p]) 
between the GSD and the CZD. This will be discussed 
further in the final section. 
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LARGE SIZE SCALING OF THE GSD AND CZD 



In Figures [7] and [8] we present the large-size behaviour 
of the GSD and CZD from the full simulations. The 
data are plotted in order to test the common large-size 
functional form suggested by the fragmentation equation 
approach for the GSD and by the GWS for the CZD , 
namely exp(— cz^) (see Eqns. ([5]) and ©). In all cases, 
the data do conform well to this functional form. In 
addition, we perfom fits to find the gradients p on these 
plots. In order to provide an estimate of the error in these 
fits, we adopt a similar strategy to that used above for 
the small-size scaling and bin the data using binwidths of 
size O.Olfc with fc = 1, 2, 20. The results of this fitting 
procedure are presented in Tables ITVl and IVl for the GSD 
and CZD respectively. 
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FIG. 8. Large-size CZD in logarithmic scale for i = 0, 1, 2 
and 3. The dashed line is the straight line fit to data. 

Once again we compare the exponents p from the 
Monte Carlo simulation data with the theoretical pre- 
dictions. For the GSD, the fragmentation equation ap- 
proach predicts values of 2q; -I- 1 for p. For the CZD, 
the fragmentation equation prediction is p = 3 for i = 
(Eqn. ([8])) and we conjecture that the values for i > 
will match those of the GSD. In contrast, the GWS pre- 
diction for the CZD is the universal value p — 2. The 
values from these theories are displayed in Tables HV] and 
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FIG. 7. Large-size GSD in logarithmic scale for i = 0, 
and 3. The dashed line is the straight line fit to data. 



1, 2 



o a = i+ 1 

" e = 20% 

d 6» = 100% 

TABLE IV. Average exponents for the large size scaling of 
the GSD using different bin-widths at coverage 6 = 20% and 
100% 



In Table IIVI we see that the fragmentation equation 
approach provides a useful point of reference to the ob- 
served large-size scaling exponents of the GSD. Again we 
see values that are bracketed by the two possible nucle- 
ation mechanisms for i = 1,2, whilst the behaviour for 
i = is a little below the predicted exponent oi p =^ 3. 
For j = 3 the data's exponent is below even that of the 
deposition-induced nucleation case. However, we have 
shown in Section IIIII above that the monomer density 
profile does not reach its saturation value in larger gaps. 
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3.108 ±0.012 


3.043 ± 0.043 


1 




2 


3.721 ± 0.020 


3.826 ±0.021 
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4.946 ± 0.029 


5.536 ± 0.033 
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5.464 ± 0.041 


6.530 ± 0.042 



" A = 2i + 3 
^ e = 20% 

e = 100% 

TABLE V. Average exponents for the large size scaling of 
the CZD using different bin-widths at coverage — 20% and 
100% 



so that the nucleation rate in these gaps is lower than 
predicted by the theory. This seems to provide a ratio- 
nal explanation for the discrepancies. 

The results in Table |V] for the large-size scaling be- 
haviour of the CZD are rather informative. We firstly 
observe that the Monte Carlo data exponents do indeed 
mirror those of the GSD in Table IIVI quite well. This 
means that the universal GWS prediction for p = 2 is al- 
ways wrong. We also see that the concrete prediction for 
i = from the fragmentation equations, namely p — 3, 
is well supported by the simulation data. 



V. SUMMARY AND CONCLUSIONS 

We have investigated one-dimensional (1-D) point is- 
land nucleation and growth simulations in order to test 
predictions for the asymptotics of the gap and capture 
zone size distributions (GSD and CZD respectively). The 
work shows that the fragmentation equation approach 
provides a good framework in which to understand the 
Monte Carlo simulation results. The theory can be used 
to investigate two cases for the nucleation process for 
i > 0, the first where nucleation is driven by deposition 
events, the second where fluctuations caused solely by 
monomer diffusion induce nucleation. 

Firstly we presented single gap simulation results 
which show that both these nucleation processes are ac- 
tive, so that the observed nucleation rates are bracketed 
by these two extremes. Furthermore, we showed that 
for larger gaps the average monomer density profile does 
not reach the long-time steady state assumed in the frag- 
mentation equations. As a result, the nucleation rates in 
large gaps are slower than predicted by the theory, with 
the shortfall increasing with gap size. Therefore, the sim- 
ple power-law scaling of the nucleation rate with gap size 
breaks down at larger sizes, with obvious consequences 
for the fragmentation equation predictions for the GSD. 

We note here that deviations from the original Black- 
man and Mulherani^ predictions for the nucleation rate 
dependence on gap size have recently been observed for 
the i = 1 1-D modelSi. In this work, the authors report 
that the nucleation rate has two regimes; for small sizes, 
it approximately obeys s*, whilst at larger sizes it approx- 



imately follows s^. The latter power-law feeds into the 
asymptotic form of the GSD and hence the CZD, yielding 
the functional form exp(— s'^). We note here that these 
values are close to those we find for i = 1 in TablcUfor the 
small gap nucleation rates and Tables IIVI and |V] for the 
large-size GSD and CZD scaling. We therefore propose 
that the explanations presented here in terms of com- 
peting nucleation mechanisms and unsaturated monomer 
density profiles will also explain the results reported in 

HI- 

We also presented data for the full island nucleation 
and growth simulation. For the small-size GSD scaling, 
we found results consistent with the fragmentation equa- 
tion predictions for i — 0. For i — 1,2,3 the exponents 
were bracketed by the values for the alternative nucle- 
ation mechanisms as expected. For the large gap size 
scaling, the Monte Carlo data followed the functional 
form suggested by the fragmentation theory, with the ex- 
ponents again being largely bracketed by the predicted 
values, although the breakdown of the nucleation rate 
scaling is apparent, especially for larger i. 

In the case of the CZD, we once again successfully 
placed the observed simulation data into the context pro- 
vided by the fragmentation equations. Interestingly, the 
GWS predictions for the small-size CZD scaling work 
extremely well since they bisect the exponents from 
the alternative nucleation mechanisms. As discussed 
elsewhere^, the predicted formula for the parameter /3 
of the GWS can be brought into line with either nucle- 
ation mechanism following the arguments of Pimpinelli 
and Einstein^^, but the original prediction of these au- 
thors (Eqn. ([TU)) ) does seem to speak well for their phys- 
ical intuition^'^. 

However, the predicted GWS form for the large-size 
CZD scaling fails badly when confronted with our 1-D 
point island simulation results. This is in contrast to re- 
cent tests performed using two-dimensional substratesii, 
which suggests that there is something unique to the 1-D 
case, possibly due to the topological constraints in how 
capture zones are constructed from the inter-island gaps. 
This aspect is worthy of further investigation. 

In order to predict the asymptotics of the CZD, we 
have assumed that the capture zones can be constructed 
from pairs of gaps sampled randomly for the GSD (see 
Eqn. (in])). This is valid provided that the nucleation has 
effectively mixed up the gaps so that nearest neighbours 
are no longer correlated^-. One consequence is that the 
small-size exponents of the CZD (say pi) are related to 
those of the GSD (say P2) through pi — 2p2 + 1. Looking 
at the results in Tables |IT] and IIIIl we see that this rela- 
tionship is reasonably obeyed for z = 0, 1 but starts to 
break down for i = 2,3. This is perhaps understandable, 
since for the higher critical island sizes, the nucleation 
rate slows down dramatically over time suggesting less 
well-mixed systems. This is another point for further 
consideration in future theory development work. 

Despite the limitations of the fragmentation equation 
approach used in this work, such as its failure to capture 
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the time-dependent nature of the monomer density pro- 
file within gaps, it has provided an excellent theoretical 
framework from which to consider the island nucleation 
process. Hence, alongside the points discussed above, 
future work might also look at how the fragmentation 
kernels can incorporate this time dependency, and how 
the two nucleation mechanisms can be combined into a 
consistent set of fragmentation equations. 
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